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Abstract 

The behaviour of neurons under the influence of periodic external input has been modelled 
very successfully by circle maps. The aim of this note is to extend certain aspects of this analysis to 
a much more general class of forcing processes. We apply results on the flbred rotation number of 
randomly forced circle maps to show the uniqueness of the asymptotic firing frequency of ergodically 
forced pacemaker neurons. The details of the analysis are carried out for the forced leaky integrate- 
and-fire model, but the results should also remain valid for a large class of further models. 

1 Introduction 

Already in 1907, long before the molecular mechanisms of neural signal transduction had been clari- 
fied, Louis Lapicque proposed a simple model for the firing behaviour of a neuron [Tl[5][3]. A crucial 
feature of this so-called integrate-and-fire model (IFM) is the separation of time-scales: the stereo- 
typical and extremely fast generation of an action potential is thought of as being concentrated in a 
single moment of time, whereas the much slower evolution of the membrane potential in the inter- 
spike intervals is modelled as a continuous process. For many questions concerning the behaviour 
of neural systems this level of abstraction turned out to be exactly the adequate one, such that 
even nowadays, more than a hundred years after Lapicque's original paper, the different variations 
of the IFM still play a central role in theoretical neuroscience 3]- One of their great achievements 
was the explanation of so-called 'paradoxical segments' that were discovered in the experimental 
investigation of pacemaker neurons in the nervous system of crayfish {Procambarus clarkii), sea 
slugs {Aplysia californica) and horseshoe crabs {Limulus polyphemus) [5l|6]. The counter-intuitive 
observation that was made in these experiments was that an increase in the frequency of periodic 
inhibitory presynaptic input can lead to an increase of the post-synaptic firing frequency. This para- 
doxon was explained by relating the respective theoretical models to monotone circle maps whose 
rotation number equals the ratio between the input and the output frequency. When the circle 
map has a stable periodic orbit, then input and output frequency remain directly proportional on 
a small neighbourhood (the paradoxical segment), disrespective of whether the input is inhibitory 
or excitatory O [71 El E] (see also Section [2]). Similar ideas have also been pioneered before by V. 
Arnold in the study of cardiac cells |10l 111] . 
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Figure 1.1: Three species for which mode-locking phenomena in the nervous system have been 
investigated experimentally From left to right: Procambarus clarkii, Aplysia californica and 

Limulus polyphemus [12) . 
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It must be said, however, that despite the great success of IFMs and the long history of their 
investigation their rigorous mathematical description is still restricted to a few very special situa- 
tions. In particular, the only types of external forcing that can be treated analytically so far are 
either periodic [13lll4j or stationary stochastic input |15II16I[T7] . Even the superposition of the two 
- noisy periodic input - is mostly accessible only by numerical methods pS]. Our goal here is take 
up the ideas used in the analysis of the periodically forced IFM an to extend these to more general 
forcing processes. Thereby, we restrict ourselves to deterministic and/or random forcing, although 
it should be possible to adapt the approach to models generated by stochastic differential equations 
as well. In order to state the main results, we first recall the construction of the IFM. 

The membrane potential V(t) of a neuron TVi remains between a lower threshold Vi and an 
upper threshold Vu- V{t) can never drop below Vi due to physiological constraints, whereas when 
it reaches Vu the neuron 'fires', meaning that an action potential is triggered and V{t) drops back 
to a rest potential Vr £ [Vj, Vu). Between the two thresholds, the potential evolves according to an 
infinitesimal law 

(1.1) V{t) = F{t,V{t)) 

with right side F : ^ R that should satisfy F{t, V;) > Vt £ R. The dependence of F on J 
corresponds to the influence of external time-dependent factors. The reset procedure when V{t) 
reaches Vu is usually expressed as 

(1.2) V{t+) = Vr if V{t) = Vu , 

where denotes the right-hand limit. Identifying the interval [V!,Kj) with the circle — R/Z, 
this gives rise to a non-autonomous circle flow (see Figure [T!T|) . 

For fixed initial values V{to) = xq, we denote by t„ the time of the n-th firing of the neuron 
A/i. Then a very basic and fundamental question is that of the existence and uniqueness of the 
asymptotic firing frequency: under what assumptions does the limit 

(1.3) i^AAi = lim n/t„ 
exist and when is it independent of the initial values to and xo? 

v(t) 
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Figure 1.2: Construction of the potential flow V{t) generated by (|l.ip and (|1.2p with initial values 
Vito) = xo. 

In the simplest case the external input is periodic in time with period p G R^. As mentioned, 
this situation is quite well-understood and has been studied for a number of different versions of 
the IFM |13l I14j . The analysis depends on the choice of a suitable Poincare section for the flow, 
by which one obtains a circle map <; : — > T'^ whose lift G : R — * R generates the rescaled 
sequence t„, that is tn+i/p — G"{tn/p). (We briefly recall the construction in Section[51) Under 
suitable assumptions on the function F this map g has good monotonicity properties that ensure 
the existence and uniqueness of the rotation number 

(1.4) p{G) = lim G"(to/p)/n = yjj\lv 

n — >oo 

and thus of the asymptotic flring frequency. As indicated above, the existence of a stable periodic 
orbit for the map g yields an explanation for the 'paradoxical segments' in situations with inhibitory 
presynaptic input. 
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A good framework to study more general forcing processes is provided by the theory of random 
dynamical systems, as exposed in |19) . In order to model the external input we will assume that 
there is an underlying forcing process, given by a measure-preserving or metric dynamical system 
(metric DS), that is, a quadruple [Q, B, /x, w) where (0, B, y,) is a probability space and : R x — > 
Q,{t,8) I— > uit{d) and tj is a flow on G that leaves the probability measure p invariant (compare 
[19]V We usually write ujt{9) — 9 ■ t. In this setting is replaced by 

(1.5) V{t) = F{eo-t,V{t)) 

where F : x R — > R and 9o is some initial value in 0. In order to ensure that (|1.5|l generates a 
RDS, we have to impose some standard technical conditions on F. We say F is uniformly Lipschitz- 
continuous in V if there exists a constant L > such that 

(1.6) \F{e,V) - F{9,V')\ < L-\V-V'\ \/9 e Q, V,V' £ R . 

For the moment, we will just assume that F is bounded and uniformly Lipschitz-continuous in V. 
These assumptions seem quite reasonable from the physiological point of view. It is also possible 
to weaken them further to some extent and we will discuss this issue in detail at the beginning of 
Section H) 

Due to the lack of a suitable structure on 0, it is not feasible anymore to take a Poincare section 
in this situation. However, it turns out that instead the flow generated by p.2p and (|1.5p can 
be analysed directly by applying a result on the fibred rotation number of randomly forced circle 
flows. (Basically, what we need is a slight modification of statements from [20; that we present in 
Section [31) This leads to the following result. 

Theorem 1.1. Suppose that the metric DS (0,S, ^,a;) is ergodic and the function F in ^1.5\) is 
bounded on O x [0, 1] and uniformly Lipschitz-continuous and non-increasing in V . Then the model 
described by \1.2}) and 11.5]) has a unique asymptotic firing frequency in the following sense: There 
exists a real number v such that for ^-almost every 9o £ <d and all xq £ [VJ, Vr) there holds 

(1.7) lim n/tn = u . 

n — >oo 

A well-known example to which this statement applies is the leaky integrate- and- fire model 
(LIFM), in which the function F takes the form 

(1.8) F{e,V) = -a{e) ■V + b(9) . 

Here a : ^ R+ corresponds to the membrane permeabilitj0 whereas 6 : ^ R refiects the 
external input. In order to apply Theorem 11.11 we have to assume that a and b are bounded. 
Shghtly weaker conditions are again discussed at the beginning of Section |3l The LIFM was first 
introduced by Stein [7] and then further investigated both theoretically and experimentally by 
Knight [S] [6]. In contrast to the so-called 'perfect integrator' or 'perfect IFM' {-^F = 0) with 
'infinite memory' used by Lapicque, it takes into account the exponential decay of the membrane 
depolarisation in time after excitation. 

The most subtle issue in the proof of Theorem 1 1.1 1 is the fact that when F takes negative values, 
then this might lead to discontinuities in the flow V{t) and to a lack of monotonicity with respect 
to the initial condition xo (see Figure [TT]) . At this point, the monotonicity of F in V is needed 
in order to recover a certain 'almost-monotonicity' property. The discontinuity does not present a 
problem in the setting of RDS, but impedes drawing further conclusions in the situation where the 
base flow tj is a uniquely ergodic0 system on a compact metric space 0. All these complications 
can be avoided by assuming that F is non-negative. In this case, we have the following. 

Theorem 1.2. Suppose that the metric DS {O, B, is ergodic and the function F in |i.5|) is 

bounded on O x [0,1], uniformly Lipschitz-continuous in V and F(6,Vr) > V& G 0. Then the 
following holds. 

(a) The model described by 1 1.2(1 and hi. 5]) has a unique asymptotic firing frequency: There exists 
a real number v such that for ^-almost every 9o £ Q and all xo £ [V;, Vu) there holds 

(1.9) lim tn/n = i/ . 

n — >oo 

^Usually the membrane permeability is chosen to be fixed, but as we will see it presents no additional 
cost to assume that it is time-dependent as well 

map or flow on a compact metric space is called uniquely ergodic, if it has exactly one invariant 
probability measure. 
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(b) If in addition Q is a compact metric space, F is continuous, F{9,Vu) > V6' G and uj 
is uniquely ergodic, then il.9\) holds for all initial conditions {8o,xo) G 6 x [V;,Ki) and the 
convergence is uniform in {9q,xo). 

This statement can, for instance, be applied to the so-called quadratic IFM, which is given by 
(fOj) and (fT3)) with 

(1.10) F{9,V) = v'' + i{e) , 

under the additional assumption that inf^ge ^(^) > va&yi{—V^,—V^}. (See [2] for the analysis 
of this model in the periodically forced case and further references). The standard example of 
a uniquely ergodic base flow would be the Kronecker flow 6 ■ t = 6 + t{uji, . . . ,LUd) with d ratio- 
nally independent frequencies uj\, . . . ,ujd, corresponding to the excitation of A^i by d independent 
pacemaker neurons. Although we will not discuss the topic in detail, we want to mention that 
mode-locking phenomena may appear in this setting whenever the asymptotic firing frequency is 
rationally related to the driving frequencies {u-i, . . . ,u)d) 

Acknowledgments. The author would like to thank K. Pakdaman for interesting discussions 
that initiated this research. This work was supported by a research fellowship (Ja 1721/1-1) of the 
German Research Council (DFG). 

2 Periodic forcing revisited 

In this section we briefly recall the analysis of the periodically forced IFM and in particular the 
construction of the circle map g mentioned in the introduction. This will also allow to discuss 
the modifications needed to extend the approach to randomly forced models. In order to simplify 
notation, we will from now on always assume that VJ = and Vu = ^■ 

Suppose that the function F in (|l.ip is p-periodic in the second variable, that is 

(2.1) F{t,V+p) = F{t,V) . 

Equivalently, we may assume that O = and ujt{9) = 6 + t/p mod 1 in (|1.5p . For simplicity, we 
also suppose that 

(2.2) F{t, Vr) > Vt G R . 

In this case, we may assume without loss of generality that Vr = Vi = 0, since the interval (V^i, V,-) is 
not accessible from outside and does not play a role in the description of the dynamics. The ODE 
(|1.1|) generates a flow 

(2.3) R^->R, {t,V)^$t{V) 

in the sense that x{t) = ^t-tgixo) is the solution of (|l.ip with initial conditions x{to) — xq. The 
first firing time ti can then be expressed as ti = to + inf{t > | "I>t-to(a;o) ~ 1}. Fixing the initial 
value xo = Vr, this allows to define a map G? : R ^ R, to ti that generates the sequence of 
spiking times tn, that is, tn+i = G{tn) Vn £ No. Rescaling G, we let G : ^ T^, x t-^ G{pt)/p. 
The periodicity assumption (|2.ip implies that G{t + p) — G{t) -f p and therefore G(t-|- 1) = G{t) + 1. 
Consequently G is the lift of a circle map <? : T^. Furthermore, using (|2.2p it is possible to 

show the map G is continuous and strictly monotonically increasing, such that g is an orientation- 
preserving circle homeomorphism. Such maps have a well-defined rotation number defined by the 
limit 

(2.4) P(G) = lim (G"(t) -t)/n = Yua {pt) / p ~ pt) / pn 

n — >oc n — >oc 

which always exists and is independent of t (see, for example, [H]). Of course, this immediately 
implies the existence and uniqueness of the asymptotic firing frequency Vfj--^ = i ■ p{G)~^ . 

There is also an alternative way of constructing the map g that yields some additional insight 
and which we want to discuss on an informal level. Suppose we let 

(2.5) F:R^^R\ {y,V)^ {l/p,F{py,V -n)) ilV£[n,n+l). 

Then F defines a vector fleld on R'^ that is invariant under integer translations and hence projects 
to a vector field on the two-torus T'^. Consequently, the flow 5' : R x R^ — > R^ generated by 
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Y{t) = F{Y{t)), Y G R^ projects to a flow * on It is now easy to see that the map g defined 
above is simply the return map of ^ to the Poincare section x {0}. (Note that (|2.2p ensures that 
this section is transversal to the flow direction.) However, from this point of view we also see that 
taking this particular Poincare section has something arbitrary. For example, one might just as well 
have chosen to define another circle homeomorphism g with lift G : R ^ R by taking the Poincare 
map of the vertical section {0} x . This simply corresponds to stopping the flow at times i„ = np. 
The asymptotic firing frequency could thus be obtained as i/j^/^ — p{G) /p. One can even arrive at the 
same conclusion without resort to Poincare sections at all. As F is constant in the first coordinate, 
the flow has skew product structure over the simple base flow ujtito) = to + t/p mod 1. In other 
words, is a periodically forced circle flow with lift Vl'. As circle homeomorphisms, such flows have 
a well-defined vertical rotation number and one obtains Ufj--^ = p{'i!)/p — limj^oo 7r2 o $t(io,a:o). 

Now, in the case of periodic forcing all this is surely a mere tautology. However, things become 
quite different as soon as one wants to consider more general types of forcing. When the driving 
space O is an arbitrary measurable space, then taking a vertical Poincare section does not make 
sense anymore, whereas the Poincare return map to O x {0} just yields a self-map of Q that might 
be very difficult to analyse. In this context, the fact that we can also directly consider the forced 
circle flow ^ generated by (|1.2p and (|1.5p turns out to be a great advantage. In the next section, 
we will present a result on randomly forced circle flows and their lifts that ensures the existence 
and uniqueness of the vertical or fibred rotation number of such skew product flows. The remaining 
sections are then dedicated to the construction of suitable lifts for the potential flow in the situation 
of Theorems 11.11 and 11.21 which immediately leads to proofs of the respective statements. 

Finally, before turning to the rigorous analysis in the next sections, we briefly want to indicate 
why further technical problems appear when condition (|2.2p does not hold anymore. First of all, this 
evidently leads to discontinuities of the flow when the potential 'jumps' from Ki = 1 to 1/r £ (0, 1) 
or, when considering lifts, from Vu + n to Vr + n+ 1 for some n G Z. This would in itself not present 
a serious problem, since most results on the existence and uniqueness of rotation numbers only 
require monotonicity, whereas continuity is less important. However, it turns out that dropping 
condition (12. 2p may at the same time lead to a lack of monotonicity of the circle flow. An example 
of how this can happen is indicated in Figure . This remains true even when F is non-increasing 
in V, as assumed in Theorem 11.11 However, in this situation it is nevertheless possible to make 
the argument work. In the case of periodic forcing, one can show that lift G of the circle map g 
constructed above is monotonically increasing on the image of one of its iterates (as done in [131114] ). 
On the level of forced flows, one may similarly show that while two different orbits may reverse 
their order, they will always remain within a bounded distance of each other. (This is the point of 
view we will adopt in Section[5l) In both cases, this is still sufficient to guarantee the existence and 
uniqueness of the rotation number. 

3 Uniqueness of the fibred rotation number 

As mentioned before, we say (0,S, is a metric DS if {Q,B,n) is a probability space and 

w : K X O O, {t,6) (jJt{9) is a flow that preserves /i, meaning that /i o cjl^^ = /i Vt G R. As 
before, we write uit{d) — 9 ■ t. A random dynamical system (RDS) over tj is a measurable mapping 

(3.1) *: RxOxR^OxR, {t,e,x) ^ {6 ■ t,%l}t(e,x)) 
that satisfies the cocycle property 

(3.2) ■ii)t^+t,_{e,x) = ipt^{e ■ti,^tAs,x)) v(ii,i2,6',x) G R^ X e X R . 

We say is a continuous RDS if for all ^ G O the mapping [t, x) h-^ ^tiP, x) is continuous. 

Evidently, what we are interested in is the existence of uniqueness of the fibred rotation number 

(3.3) /9(VE') = Vixa{%l}t{e,x)-x)/t 

t — ^oo 

when is obtained as the lift of a randomly forced circle flow. For this, we will use the following 
assumptions: 

• There holds 

(3.4) iljt{0,x + l) = iljt{0,x) + l V(t,6',3;) G T X e X R . 

In particular, this implies that projects to a random dynamical system on the circle T^. 
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Figure 2.1: Lack of monotonicity of the lift of tiie potential flow. Two solutions Vi{t) and V2{t) 
with Vi{to) — xi < X2 = V2(to) have changed order at time s. 



• There exists a constant C > such that 

(3.5) M^,^) < V't(f,a;') +C W{t,6,x,x') eTT^e xR^ with X < x' . 

In other words, the mappings ipt{0, .) are 'almost monotone' in x, up to a uniform constant C. 

• There exists a constant > such that 

(3.6) \i^t{e,x)~x\ < ri y{t,6,x) £[0,1] X e xR . 

Theorem 3.1. Suppose that {0,B, ij,,uj) is ergodic and the RDS ^ satisfies ^3. 4^ - 113. 6\) . Then there 
extsts a real number p and a set Oo C B of full measure, such that 

(3.7) \im{ip{t,e,x)~x)/t = p y{e,x) eOo xR . 

t — 'OO 

Remark 3.2. (a) When is a compact metric space and ui a uniquely ergodic flow on G, then 
this result is well-known and due to Herman [23]. (See also [2^ for a precursor in the context of 
quasiperiodic Schrodinger operators.) For the more general case of ergodically forced monotone 
circle maps, it was proved more recently by Li and Lu [20J . The only new aspect here is that 
we work with continuous time replace monotonicity by the slightly weaker property p.Sp . It is 
not surprising that this can be proved along the same lines as the previous results, but since the 
proof is rather short anyway we include the details for the convenience of the reader. 

(b) When is a continuous RDS on the circle, then assumption (|3.5p can be replaced by the much 
more general one that 

(3.8) lira \'ipt{e,x) -tpt{e,x')\/t ^ V(6', a;, a;') G 6 x . 

t — >oo 

In this case Herman's original proof, which uses the existence of a 'I'-invariant probability measure 
that projects down to /i, remains valid with hardly any modifications. However, since we do not 
assume 'I' to be continuous (and this is crucial for the application to the forced LIFM), such an 
invariant measure does not necessarily exist. 

Proof of Theorem\3Ji By and we have 

(3.9) \{ipt{d,x)-x)-{iPt{9,x')-x')\ < C+1 \/{t,e) gT xe,x,x' € Ag . 

(Note that due to (|3.4[) we may assume w.l.o.g. that x £ [x' — l,x')), in which case the estimate 
is a direct consequence of p.Sp .l It follows that when the limit p{9,x) := \imt-^oo{i>t{d , x) — x)/t 
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exists for one a; £ R, it exists for all a; £ R and does not depend on x. Given any n € No, let 
F^{e) ■- ipni0,O) + C + 1. Then implies that 

ED 

< i)m{e ■ n, 0) + C + 1 + Fn{e) = F„(6i ■ n) + F„{e) . 

Thus the random variables F„ form a subadditive sequence over the measure-preserving transfor- 
mation 9 9 ■ 1. Further (|3.6p implies that Fi is bounded. Hence, we can apply Kingman's 
Subadditive Ergodic Theorem (see, for example, [H]), which yields the existence of a ^-integrable 
function p : O ^ R and a set 9' C O, such that for all 9 £ Q' there holds 

lim F„{9)/n = p{9) . 

n — ^oo 

(|3.6p and (|3.9[) together now now imply that for all {9, x) G 0' x R we have 

lim {ipt{9,x) -x)/t = p{9) . 

t — *oo 

It is easy to see that the function p is invariant, that is p{9 ■ t) = p{9) \/t G T, and the ergodicity of 
the base flow uj therefore implies that p{9) — p dp, —: p on a, set of full measure O". If we now 
let Oo = O' n O" then all the assertions of the theorem are satisfied. □ 



4 Construction of the potential flow 

The aim of this section is to formalise the model described by (|1.2|l and (|1.5|1 in the introduction and 
to prove Theorem 11.21 (The proof of Theorem 1 1 . II will then be given in Section O) More precisely, 
we will construct a lift V for the circle flow V that describes the evolution of the membrane 
potential (recall that we identify the interval [V;,l/u) with the circle. We assume without loss 
of generality that VI = and Vu = 1. Unfortunately, the construction has a somewhat technical 
flavour, which basically comes from the need to treat the discontinuities produced by (|1.2p in a 
formally precise way. However, in order to treat the problem in a rigorous way a certain amount 
of detail seems unavoidable, in particular since the 'almost-monotonicity property' needed for the 
proof of Theorem 1 1.1 1 is a rather subtle issue. 

Suppose that {Q,B, p,u)) is a metric DS and consider the non-autonomous differential equation 
(4.1) x{t) = F{9-t,x{t)) 



with _F : R X S ^ R, which is equivalent to (|1.5|l . As mentioned in the introduction, we first discuss 
the precise technical conditions on F needed for the construction. Given any / : R ^ R, we let 

ll/ll[o,i],o,o = sup |/(2:)| and 

a:e[0,l] 

llfll - llfll -U «„T^ I/fa) - 

ll/ll[o,i],o,i - ll/ll[o,i],o,o + sup : . 

x^yelO.l] \y ~ ^\ 

(The terminology follows 119} Appendix B].) We then require that there exists a constant C > 0, 
such that for all S £ there holds 

(4.2) / \\F{9 ■ t,. )\\io,i],o,o dt < C and 

Jo 

(4.3) / \\F{9 ■ t,. )\\[o,i],o,i dt < oo. 

Jo 

The restriction to the interval [0, 1] in the definition of the norms is explained by the fact that due 
to the reset procedure described by (|1.2|l this is the only part of the phase space we are interested in. 
We could equally assume that the conditions are satisfied on all of R and just modify the function 
F if they are violated outside of [0,1]. Assumption (|4.2|) is needed to ensure that the hypothesis 
of Theorem 13. II in the appendix are met - roughly spoken, it just ensures that action potentials or 
spikes cannot be generated at an arbitrary rate (there will be at most C/(l — Vr) + 1 spikes in any 
interval of length 1). (|4.3|) is just a standard technical condition which ensures that (|4.ip generates 
a RDS 

(4.4) $: RxOxR^exR, ${t,9,x) ^ {9 ■ t,ipt{9,x)) 
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over the base flow uj (see [191 Theorem 2.2.2] )|!] From now on, we will always assume that $ is the 
skew product flow generated by (I4.1|l . 

In the following, we will give a precise definition of the lift of the potential semiflow 

(4.5) V : R+xexR, {t,e,v) ^ {6 ■ t,Vt{e,x)) 

corresponding to the model described by p.2p and (|1.5|l . The membrane potential at time t with 
initial values t = Q, 6 and xo is then obtained as V(i) — Vt{6,xo) mod 1. The only assumption on 
F needed for this construction is (|4.3|l . (|4.2|l and the additional assumptions made in Theorems 11.11 
and ll.2l will be required only later in order to ensure that V meets the assumptions of Theorem 13. II 
Given any a; £ R, we denote by [x] £ Z its integer part and by {x} = x—[x] G [0, 1) its fractional 
part. For any {6,x) £ O xM. we define 

(4.6) T{e,x) = mi{t>0\(pt{e,{x}) ^1} 

with t(6, a;) = oo if the set on the right is empty. This is the time when the first spike is generated. 

For any s G (R+)'* we denote by S = S{s) the sequence given by Sn ~ X^iLo Then, given 
(6,x) G X R, we recursively define a sequence s„ = Sn{9,x) and the corresponding sequence 

s„ = Sn{e) by 

(4.7) So = , si = r(6l, x) and s„+i = r(6l ■ S^, K) . 

(Recall that Vr G [0, 1) = \Vi, Vu) is the rest potential introduced in (|1.2p .'l Sn is the time when the 
n-th action potential is triggered, whereas s„ is the length of the time interval between the n — 1-th 
and the n-th spike. Given [t, 8, x) G R^ x x R we define n{t, 9, x) as the unique integer n such 
that Sn{0,x) <t< S„+i{9,x), that is 

(4.8) n{t,e,x) = max{n G N j S„(6l,x) < t} . 

In other words, n{t, 9, x) is just the number of spikes generated until time t. The lifted potential 
flow is then defined by (14. 5|) with 



(4.9) Vt{9,x) 



[x] + ipt{9,{x}) iin{t,9,x) = ■ 

[x] + n + ipt-s„ (9 ■ S„, Vr) if n{t, 9,x) =n . 



Remark 4.1. (a) The circle flow V defined by (|1.2p and (ll.5|l is typically not continuous and 
therefore a priori has many different lifts that do not only differ by an integer constant. However, 
the particular lift V defined above is the unique one that 'counts' the number n{t, 9, x) of spikes 
(action potentials) generated up to time t, in the sense that n(t, 9, x) equals the number of integers 
in the interval [x, Vt{9, x)]. In this way, the asymptotic firing frequency is obtained as the inverse 
of the rotation number 

(4.10) p{V) = lim {Vt{9,x))/t , 

t — >oc 

provided this limit exists. To show that this is the case under the assumptions made in Theo- 
rem [L2] is the aim of the remainder of this section. 



(b) We note that 

(4.11) Vti9,x + 1) = Vtie,x) + 1 V(f, 6», x) G R+ x e X 

Further, when (14.211 holds then 



(4.12) \Vt{9,x)~x\ < -^-^ + 1 Vt G [0,1], (6l,x) G e X R . 

Finally, we note that by definition we have 

(4.13) VsJ9,x) = [x]+n + Vr . 

Consequently the mapping t Vt{9,x) is continuous from the right. The only discontinuities 
occur exactly at times t = S„, except in the case K- = in which t >—f Vt{9,x) is continuous. 



^Strictly spoken, we would have to take into account finite escape times, such that t i— » ipt{9,x) is only 
defined on a subinterval of R. However, since we are only interested in the dynamics on X [0, 1] we assume 
for simplicity that 114. l| l always has global solutions. 
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Proposition 4.2. The mapping V given by and _?0| ) defines a RDS over the base flow u). 

In particular, it has the cocycle property 

(4.14) Vt,+t2{0,x) = Vt^{e-ti,Vt,{e,x)) vti,t2eR+, {e,x)eexR. 

Furthermore, if Vr = and F{6,0) > G then the mapping x ^ 'Vt{Q,x) is monotonically 
increasing for all {t,0) £ R"'". // in addition u is a continuous flow on a topological space O, F in 
114. j|] is continuous and F(6, 1) > V6 G 0, then V is continuous as a mapping xOxR— > GxR. 



This shows that in the situation of Theorem 11.21 the semi-flow V satisfies all assumptions of 
Theorem 13.11 which proves part (a) of the theorem. Part (b) for uniquely ergodic base flows then 
follows from Herman's classical result on the fibred rotation number [251 Theorem 5.4]. 

Proof. According to the definition of a RDS, we have to show that the mapping V is measurable 
and has the cocycle property. (Since measurability follows from standard arguments we omit the 
details.) 

Cocycle property. In order to show (|4.14ll we need the following technical statement. 

Claim 4.3. Given t\,t2 > and {6,x) G O x R, let ni = n{t\,9,x), a\ = Sni{p,x), n2 = 
n{t2,6 ■ ti, ((9tj_o-i {6 ■ ai,Vr)), 0-2 = S„2 {6 ■ ti, ipti^ai {d ■ ai, Vr)) and n — n{ti +t2,9,x). Then 

(4.15) n = ni + n2 and 

(4.16) S„{e,x) = ti+aa. 

Proof. First, we note that ni = n{t\,9,x) = n{ai,8,x) and n2 — n{t2,0 ■ ai,Vr). Hence, we 
may assume without loss of generality that ti = cri. We claim that 

(4.17) Si:{e-s,{9,x),Vr) = Sk+,{e,x)-Sj{e,x) vj.fcGN. 

In order to see this, we proceed by induction on k. There holds 

Si{e-Sj{6,x),Vr) = T{6-Sj{e,x),Vr) = sj+i{e,x) = Sj+i{e,x) ~ Sj{9,x) . 

Further, if (|4.17p holds for k, then 

Sk + l{S ■ Sj{9,x),Vr) = t{0 ■ {Sj{e,x)+Sk{0 ■ Sj{6,x),Vr)),Vr) + Sk{0 ■ Sj{6,x),Vr) 

= T{e- Sk+j {9,x),Vr) + Sk+j {e,x)-Sj{0, x) 

= sk+,+i{e,x) + Sk+,{o,x)-s,{e,x) = Sk++j+i{e,x)-Sj{e,x) . 

This proves (|4.17|l . Now, we have 

n{t2,e ■ ai,Vr) = ma.x{k en \ Sk{0 ■ c7i,Vr) <t2} 

max{fc G N | Sn-i^+k{d,x) — ai < t2} 
= max{fc G N I S'„j+fe(6',a;) < t2 + o-i} 
= m£tx{Sk{d,x) < t2 + ai} ~ ni = n{t2 + o-i,6,x) — ni . 

Thus l|4.15p holds. In order to show (|4.16|l . recall that we assumed ti — ai. Therefore 
ti + aa = ai+S„,ie-ai,Vr) S„{e,x) . 

o 

Now the cocycle property follows easily. Given ti,t2 > 0, we define ni,n2,(Ti and 02 as in the 
claim above. First, suppose that ni,n2 > 1. Then 

Vt,{e-ti,Vt,{e,x)) [x] + ni + Vt,{e-ti,vt^-ad9-'^i,Vr)) 

= [x] + n-i + n2 + ipt2-<T2{9 ■ {ti + '^2),Vr) 

[x] + n + ipt^+t2-s„(e,x)iS ■ Sn{0,x),Vr) = Vti+t2{6,x) . 
The cases where ni = or n2 = are treated similarly. 
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Continuity and monotonicity. Suppose that F is strictly positive and Vr = 0. As mentioned 
in Remark 14.11 the mapping {t,x) i-^ Vt{8,x) is continuous in t in this case. In order to show the 
monotonicity, fix 61 G 9 and xi < X2 eM. and let to = mi{t G R+ | Vt{9, Xi) > Vt{9, X2)}. 

Suppose for a contradiction that to < 00. Then continuity in t yields Vq := Vto(^ia^i) = 
Vt„{9,X2)- We distinguish two cases. First, assume that Vo ^ Z + Vr- In this case the orbits of 
{6,xi) and {6,X2) coincide with integer translates of orbits of the flow "1> on a small interval Jo 
around to- They are therefore either equal or distinct on all of Jo, but cannot merge exactly at time 
to. Secondly, assume that Vo ~ k + Vr for some k £ Z. Since F is strictly positive, this would imply 
that limt^tp Vt{8,xi) — limt^to Vt{8,X2) = k, which would again mean that two distinct orbits of 
the flow "1> would have to merge at time to- Hence, in both cases we arrive at a contradiction. 

Now assume in addition that a; is a continuous flow on a topological space O, F is continuous 
and 

(4.18) F(e,i) > V6» G e . 

Then F is bounded on G x [0, 1], which together with the uniform Lipschitz continuity oi F in V 
implies that for any compact interval ICR the flow "I> generated by (|4.ip is uniformly continuous on 
7 X O X [0, 1]. Now, the flow V is obtained by concatenating integer translates of finite trajectories of 
<1> in X [0, 1]. V will therefore inherit the uniform continuity of "1>, provided that no discontinuities 
are created by this concatenation in (|4.9p . 

In order to see this, observe that (I4.18P together with the continuity of F implies that r defined 
in (|4.6p is continuous in {9,x), except when a; G Z. By induction, this yields that the spiking times 
S„{9, x) depend continuously on {9, x) as well unless a; G Z. Consequently n{t, 9, x) is continuous in 
{t,9,x) unless x £ Z or t — Sn{6,x) for some n G Z. Furthermore [x] is obviously locally constant 
when X ^ Z. Thus, it only remains to check that Vt{9,x) defined by (|4.9p is continuous in {t,9,x) 
when a; G Z or when t = Sn{9,x). However, this can be seen quite easily by having a careful look 
at (|4.9p . When x passes an integer, then [x] will jump up by one, but at the same time n{t,9,x) 
will drop down by one, such that the two discontinuities cancel each other. Similarly, when t and 
Sn{9, x) change order then n{t, 9, x) has a discontinuity of size 1, but at the same time the quantity 
ft-s^ft e ^)[9 ■ S„(t,0,x)[9,x),Vr) jumps by one in the opposite direction. □ 

5 Modifications needed for the forced LIFM 

Proposition 5.1. Suppose that the function F in j| ) is non-increasing in x. Then the mapping 
V defined by ^.5^ ) and satisfies 

(5.1) Vt{9,xi) < Vt{9,X2) + K \f{t,9,xi,X2) £9.^ X e xR"^ with xi < X2 , 

where K — + 1. 

Proof We first prove that for all (t, 0) G R+ x 9 there holds 

(5.2) Vt(9, X2)-l < Vt{9,xi) < Vt{9,X2) Va;i, 2:2 G [K, 1) with xi < 2:2 . 

Fix 61 G 9 and xi < X2 £ [Vr,l). Let S„ := 5'„(6l,xi) and S'„ := Sn(9,X2). Then we show by 
induction on n that 

(5.3) Vt(9,X2)-l < Vt{9,xi) < Vt{9,X2) yt£[0,S'„] and 

(5.4) Vs',S^,xi) £ (n-l + K-,n). 

In order to do so, we use the fact that for all 9 £ O, t > and x < x' there holds 

(5.5) <pt{9,x') — ipt{9,x) < X —X 

since ^ 

H^t(9,x)-ipt(9,x) = X -x^ [ F{9 ■ r,ipr{9,x')) - F{9 ■ r,ipr{9,x)) dr 

Jo 

and the integral on the right is non-positive since F is non-increasing. In order to start the induction, 
note that t{9,xi) > t{9,X2). Hence, if t < 5*1 = t{9,X2), then 



(5.6) 



< Vt{9,xi) =ifit{9,xi) < ift{9,X2) 



= Vt{9,X2) < 1 
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Further (|5.5ll . a;2 — a;i < 1 — 14 and the fact that lim^^g/ Vt{0, X2) ~ 1 imply that 

Vr < Vs;(e,Si) < 1 . 

This shows that ([0|) and (fO)! hold for n = 1. 

Now, suppose that (|5.3|l and H5.4|) hold for n > 1. Then using (|5.5p and limt/s,^ Vt(S,a;i) = n 
in order to compare Vt{6,X2) — 1 and Vt{0,xi), similar as in (15. 6p . we obtain that 

n-1 < Vt{e,X2)-l < Vt{e,xi) < n VteK.Sn) and 
n-l + Vr < Vs„{e,X2) - 1 < n . 

This shows that (|5.2|l holds on the interval [S'n,Sn], and for the remaining interval {Sn,S'n+i] we 
can now proceed in exactly the same way as in the case n = 0. This proves (|5.3p and (|5.4p for all 
n G N and hence (|5.2p . 

In order to treat the general case, now assume that xi < X2 are arbitrary. Due to (14. lip we may 
assume without loss of generality that xi £ [0, 1) and X2 € {xi,xi + 1). We distinguish two cases. 
First, assume that X2 G {xi, 1). Let S'„ = Sn{6,X2) as above and 

m ■- max{n e N | Vs'^i9,xi) G (0,14)} • 

Then, since Vs'^{9,xi) — ips'^{fi;X\) is still below Vr ~ Vs!^{0, X2) — m we have that Vg/ ^^{9,xi) 
is still below 1. Consequently both Vs' ^{d,xi) and Vg' ^{0,X2) are contained in [14,1). Using 
the cocycle property l|4.14p we can therefore apply (|5.2p to see that 

Vtie,xr)-i < Vt{e,x2)-{m + i) < Vt{e,xi) •^t>sL+i- 

We thus obtain 

(5.7) Vt{9,X2)- (m + l) < Vt{e,xi) <Vt{e,X2) Vt G R+ . 

This already proves (15. ip for such a;i,a;2. 

In order to treat the case where X2 G (1, 1 + 3:^1) we have to obtain some information on m. More 
precisely, we claim that 

(5.8) ^<_^ + l. 
This follows from the fact that for all n < m there holds 

Vs'je,xi) > (n-i).(i-v;) 

which can be proved easily by induction using (|5.5p together with the fact that Vg' {9, X2) ~n = Vr 

and lim^ ygi Vt{0, X2) — n = 1. 

Now assume that X2 G (1, 1 + xi). Then we can apply (|5.7p to a;i = a::2 — 1 and x'2 = xi to 
obtain that 

Vt{9,xi)~{m + 1) < Vt{9,X2)-l VtGR+ 

and hence 

Vt{9,xi) < Vt{9,X2)+m VtGR+. 
Together with (|5.8p this proves (15. ip . □ 
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